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<C '. 

We study smoothness of densities for the solutions of SDEs whose 
• coefficients are smooth and nondegenerate only on an open domain D. 

We prove that a smooth density exists on D and give upper bounds for 
this density. Under some additional conditions (mainly dealing with 
^■y^ the growth of the coefficients and their derivatives), we formulate up- 



per bounds that are suitable to obtain asymptotic estimates of the 
density for large values of the state variable ("tail" estimates). These 
' results specify and extend some results by Kusuoka and Stroock 

[J. Fac. Sci. Univ. Tokyo Sect. IA Math. 32 (1985) 1-76], but our ap- 
proach is substantially different and based on a technique to estimate 
the Fourier transform inspired from Fournier [Electron. J. Probab. 
13 (2008) 135-156] and Bally [Integration by parts formula for lo- 
cally smooth laws and applications to equations with jumps I (2007) 
^ ' The Royal Swedish Academy of Sciences]. This study is motivated 

oo : by existing models for financial securities which rely on SDEs with 

in . 

non-Lipschitz coefficients. Indeed, we apply our results to a square 
root-type diffusion (CIR or CEV) with coefficients depending on the 
state variable, that is, a situation where standard techniques for den- 
00 ' sity estimation based on Malliavin calculus do not apply. We establish 

, the existence of a smooth density, for which we give exponential es- 

timates and study the behavior at the origin (the singular point). 



1. Introduction. It is well known that Malliavin calculus is a tool which 
allows, among other, to prove that the law of a diffusion process admits 
a smooth density. More precisely, if one assumes that the coefficients of an 
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SDE are bounded C°° functions with bounded derivatives of any order and 
that, on the other hand, the Hormander condition holds, then the solution 
of the equation is a smooth functional in Malliavin's sense, and it is non- 
degenerate at any fixed positive time. Then the general criterion given by 
Malliavin [15] allows one to say that the law of such a random variable is 
absolutely continuous with respect to the Lebesgue measure, and its density 
is a smooth function (see [16] for a general presentation of this topic). 

The aim of this paper is to relax the aforementioned conditions on the 
coefficients: roughly speaking, we assume that the coefficients are smooth 
only on an open domain D and have bounded partial derivatives therein. 
Moreover, we assume that the nondegeneracy condition on the diffusion 
coefficient holds true on D only. Under these assumptions, we prove that 
the law of a strong solution to the equation admits a smooth density on D 
(Theorem 2.1). Furthermore, when D is the complementary of a compact 
ball and the coefficients satisfy some additional assumptions on D (mainly 
dealing with their growth and the one of derivatives), we give upper bounds 
for the density for large values of the state variable (Theorem 2.2). We 
will occasionally refer to these aymptotic estimates of the density as "tail 
estimates" or estimates on the density's "tails." 

Local results have already been obtained by Kusuoka and Stroock in [13], 
Section 4. Here the authors work under local regularity and nondegeneracy 
hypotheses too, but the bounds they provide on the density are mostly signif- 
icant on the diagonal (i.e., close to starting point) and in the small time limit, 
while they are not appropriate for tail estimates. Moreover, the constants 
appearing in the estimates are not explicit (cf. (4.7)— (4.9) in Theorem 4.5 
and the corresponding estimates in Corollary 4.10, [13]). In the present pa- 
per, we provide upper bounds that are suitable for tail estimates, and we 
find out the explicit dependence of the bounding constants with respect to 
the coefficients of the SDE and their derivatives. Our bounds turn out to be 
applicable to the case of diffusions with tails stronger than gaussian. This is 
the case for square-root diffusions, which are our major example of interest 
(see Section 3). Also, our approach is substantially different from the one 
in [13]. In particular, we rely on a Fourier transform argument, employing 
a technique to estimate the Fourier transform of the process inspired from 
the work of Fournier in [8] and of Bally in [2] and relying on specifically- 
designed Malliavin calculus techniques. We estimate the density pt(y) of the 
diffusion at a point y € D performing an integration by parts that involves 
the contribution of the Brownian noise only on an arbitrarily small time 
interval [t — 6, t) . This allows us to gain a free parameter 5 that we can even- 
tually optimize, and the appropriate choice of 5 proves to be a key point in 
our argument. We do not study here the regularity with respect to initial 
condition (which may be the subject of future work). 
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Our study is motivated by applications to Finance, in particular by the 
study of models for financial securities which rely on SDEs with non-Lipschitz 
coefficients. As it is well known, a celebrated process with square-root diffu- 
sion coefficient was proposed by Cox, Ingersoll and Ross in [5] as a model for 
short interest rates and was later employed by Heston in [10] to model the 
stochastic volatilities of assets. The stochastic-a/3/3 or SABR model in [9] is 
based on the following mixing of local and stochastic volatility dynamics: 



where < p < 1, B t = pW} + y/l - p 2 W?, (W 1 , W 2 ) is a standard Brownian 
motion, and p G [—1, 1] is the correlation parameter. In this paper, we apply 
our results to one-dimensional SDEs of the form 



where a G [1/2, 1) and a, b and 7 are C£° functions. When the coefficients a, 
b, a are constant, the solutions to this class of equations include the classical 
CIR process (a = 1/2) and a subclass of the CEV local volatility diffusions 
(when a = and b = —r). As pointed out by Bossy and Diop in [4], SDEs 
with square-root terms and coefficients depending on the level of the state 
variable arise as well in the modeling of turbulent flows in fluid mechanics. 
It is well known that for a CIR process, the density of Xt is known explic- 
itly. The main contribution of our results lies in the fact that they apply 
to the more general framework of SDEs whose coefficients are functions of 
the state variable, thus when explicit computations are no longer possible. 
Theorem 2.2 directly applies to (1.1) and allows one to show that Xt ad- 
mits a smooth density on (0, +00): under some additional conditions on the 
coefficients (mainly dealing with their asymptotic behavior at 00 and zero), 
we give exponential-type upper bounds for the density at infinity (Proposi- 
tion 3.3) and study the explosive behavior of the density at zero (Proposi- 
tion 3.4). The explicit expression of the density for the classical CIR process 
shows that our estimates are in the good range. 

The paper is organized follows: in Section 2 we present our main results 
on SDEs with locally smooth coefficients (Section 2.1), and we collect all 
the technical elements we need to give their proofs. In particular, in Sec- 
tion 2.2 we recall the basic tools of Malliavin calculus on the Wiener space, 
which will be used in Section 2.3 to obtain some explicit estimates of the 
L 2 -norms of the weights involved in the integration by parts formula. This 
is done following some standard techniques of estimation of Sobolev norms 
and inverse moments of the determinant of the Malliavin matrix (as in [16] 
and [6], Section 4), but in our computations we explicitly pop out the depen- 
dence with respect to the coefficients of the SDE and their derivatives. This 




(1.1) 
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further allows us to obtain the explicit asymptotic estimates on the density. 
Section 2.4 is devoted to the proof of the theorems stated in Section 2.1. 
We employ the Fourier transform argument and the optimized integration 
by parts we have discussed above. Finally, in Section 3 we apply our results 
to the solutions of (1.1). 

2. Smoothness and tail estimates of densities for SDEs with locally smooth 
coefficients. 

2.1. Main results. In what follows, b and o~j are measurable functions 
from R™ into R m , j = 1, . . . , d. For y G R m and R > 0, we denote by B R (y ) 
[resp., B_ R (y )] the open (resp., closed) ball B R (y ) = {y G R m : \y - yo\ < R} 
[resp., B R (yo) = {y G R m : \y — yo\ < R}], where | • | stands for the Euclidean 
norm. We follow the usual notation denoting C£°(A) the class of infinitely 
differentiable functions on the open set A C R m which are bounded together 
with their partial derivatives of any order. For a multi- index a G {1, . . . , m} k , 
k>l, d a denotes the partial derivative -m — ^ — . 

Let < R < 1 and y G R m be given. We consider the SDE 



for a finite T > and x G R m , and assume that the following hold: 
(HI) {local smoothness) b,aj G C£°{B 5R (y );R m ); 

(H2) (local ellipticity) aa*(y) > c y0iR I m for every y G B^ R (yQ), for some < 

(H3) existence of strong solutions holds for the couple (6, a). 

Let then (X t ;t G [0, T]) denote a strong solution of (2.1). Our first main 
result follows. 

Theorem 2.1. Assume (HI), (H2) and (H3). Then for any initial con- 
dition x G R m and any < t < T, the random vector Xt admits an infinitely 
differentiable density pt, yo on B R (yo). Furthermore, for any integer k>3 
there exists a positive constant depending also on y$,R,T, m,d and on 
the coefficients of (2.1) such that, setting 




(2.1) 



te[0,T],i 



= l,...,m 



P t {y) = F(ini{\X s - y\ : s G [(* - 1) V t/2, t]} < 3R) 



then one has 



(2.2) 
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for any y G Br(j/o)- Analogously, for every a G {1, . . . ,m} k , k>l, 

(2.3) \d aP t,y (y)\ <Pt(y )(l+ tm (J + 3) /2 ) A ^ + 3 

/or every y G B R (y ). 

The functional dependence of with respect to yo, R, T and to the 
bounds on the coefficients b and a is known explicitly. We provide the ex- 
pression of Afc in Section 2.4 in a more detailed version of Theorem 2.1 
(Theorem 2.4) which we do not give here for the simplicity of notation. 

When the coefficients of (2.1) are smooth outside a compact ball and 
have polynomial growth together with their derivatives therein, according 
to Theorem 2.1 a smooth density exists outside the same compact set, and 
one can deduce some more easily-read bounds on the tails. More precisely, 
we consider the following assumptions: 

(HI') There exist r) > such that b,aj are of class_C°° on R m \ B^O), and 

(H2) holds for any R > and y such that B 3R (y ) Ct m \ B v (0); 
(H4) there exist q, q > 0, and positive constants < Cq < 1 and C k , k > 1, 
such that for any a G {1, . . . , m} k 

(2-4) \daV(y)\ + \d a <T i j (y)\<C k (l + \y\«) 
and 

(2.5) aa*(y)> ColyP/m 
hold for \y\ > rj. 

Theorem 2.2. Assume (Hi') and (H3). 

(a) For any initial condition x G M. m and for any < t < T, Xt admits 
a smooth density on M m \ B^iO). 

(b) Assume (H4) as well. Then estimates (2.2) and (2.3) hold with R=l and 

(2.6) A k = A k (y ) := C k)T (l + \yo\ q ' k(q) ), 

for every \yo\ > r/ + 5. The value of the exponent q' k {q) is explicitly known 
(and provided in Theorem 2.5). 

(c) If moreover sup 0<s<t \X s \ has finite moments of all orders, then for 
every p > and k>l there exist positive constants C kiPi T such that 



(2.7) 



lft(y)l<C 3lPir (l + ^J|y| 
\daPt(y)\<c k , p , T (i + tmi2 l +3)/2 ) 



for every <t <T and every \y\ > rj + 5. 

In the above, the C k , P ,T are positive constants depending on k,p,T and also 
on m,d and on the bounds (2.4) and (2.5) on the coefficients. 



The proofs of these results will be given in Section 2.4. 
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Notation. Through the rest of the paper, (•, •) will denote the Euclidean 
scalar product in M. m , while the notation | • | will be used both for the absolute 
value of real numbers and for the Euclidean norm in M m . Furthermore, when 
= 0\, . . . , 9 V is a family of parameters, unless differently specified by Cq, 
we denote a constant depending on the #j's but not on any of the other 
existing variables. All constants of such a type may vary from line to line, 
but always depend only on the O^s. For functions of one variable, the A;th 

derivative will be denoted by 4-^ or f^ k \ We will follow the convention of 
summation over repeated indexes, wherever present. 

2.2. Elements of Malliavin calculus. We recall hereafter some elements 
of Malliavin calculus on the Wiener space, following [16] . 

Let W = (Wf, . . . , Wf; t > 0) be a d-dimensional Brownian motion defined 
on the canonical space (O, J 7 , IP). For fixed T > 0, let T~L be the Hilbert space 
U = L 2 ([0,T};R d ). For any h G U we set W(h) = £j =1 J Q T h? \s) dWl ' , and 
consider the family S C L 2 (Q,, F, P) of smooth random variables defined by 

S = {F : F = f(W(h!), W(h n )); h u ...,h n en;fe C™ ol (R n );n > 1}, 

where C^j denotes the class of C°° functions which have polynomial growth 
together with their derivatives of any order. 

The Malliavin derivative of F £ S is the d-dimensional stochastic process 
DF = {plF, D?F; r G [0, T\) defined by 

D r F = E 7T" • • • > W(h n ))hi(r), j = l,...,d. 

i=l 

For any positive integer k, the kth. order derivative of F is obtained by iterat- 
ing the derivative operator: for any multi-index a = (a\, . . . , G {1, . . . , d} k 
and (n, . . . , r k ) G [0, T] k , we set Dg^ffiF := D% ■ ■ ■ D?£F. Given p > 1 and 
positive integer k, for every F G S we define the seminorm 



l/p 



\ k:P = E[|Fn + ^E[||L>WF||^] 



where 



1/2 



and the sum is taken over all the multi- indexes a = (a\, . . . , a k ) G {1, . . . , d} k . 
We denote with B fc ' p the completion of S with respect to the seminorms 
|| • ||fe )P , and we set D°° = f] p>1 f\>i ^ k,p - We may occasionally refer to ||fc,p 
as the stochastic Sobolev norm of F. 
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In a similar way, for any separable Hilbert space V we can define the 
analogous spaces H> k,p (V) and D°°(l A ) of V-valued random variables with 
the corresponding || • ||fc jPj y semi- norms (the smooth functionals being now 
of the form F = Y^j=\ Fj v j> where Fj E S and Vj E V). In particular, for any 
Revalued process (u s ;s < t) such that u s E B fc ' p for all s E [0,t] and 

k 

\\u\\ H + ^II^^^H-H®^ 1 < °°> P " a - S -' 
h=l 

we have 

/ k \ 1 /P 

\\u\\ k , p , H = U[\\u\\ p n ] +^E[\\D^u\\^ h+1 }\ . 

Finally, we denote by 5 the adjoint operator of D. 

One of the main applications of Malliavin calculus consists of showing 
that the law of a nondegenerate random vector F = (F , . . . , F m ) E (D°°) TO 
admits an infinitely differentiable density. The property of nondegeneracy, 
understood in the sense of the Malliavin covariance matrix, is introduced in 
the following: 

Definition 2.1. A random vector F = (F 1 , ... ,F m ) E (D°°) m , m > 1, 
is said to be nondegenerate if its Malliavin covariance matrix op, defined by 

(a F ) itj = (DF\DF^n, i,j = l,...,m, 

is invertible a.s. and moreover 

E[det(a F )- p ] < oo 

for all p > 1 . 

The key tool to prove smoothness of the density for a nondegenerate 
random vector is the following integration by parts formula (cf. [16]). 

Proposition 2.1. Let F = (F 1 ,...,F m ) E (D°°) m ; m>l, be a nonde- 
generate random vector. Let G E D°° and <fi E C™](]R m ). Then for any k>l 

and any multi-index a = (ax, . . . ,ctk) E {1, . . . ,m} k there exists a random 
variable H a (F,G) E D°° such that 

(2.8) E[d a <l>(F)G\=E[<t>(F)H a (F,G)], 

where the H a (F,G) are recursively defined by 

H a (F,G) = H( ak )(F, H( au ... t a k _ 1 )(F, G)), 

in. 

H {l) {F,G) = Y,5(G{a- F 1 ) h3 DFi). 
i=i 
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2.3. Explicit bounds on integration by parts formula for diffusion pro- 
cesses. The notation of this section is somehow cumbersome, as we try 
to keep our bounds as general and as accurate as possible. The framework 
will nevertheless considerably simplify in Section 2.4, when we will give the 
proofs of the results stated in Section 2.1. 

Throughout this section, X = (Xt]t > 0) will denote the unique strong 
solution of the SDE 

rt d rt 

(2.9) X\ = x l + / B\X S ) ds + y^ A *j( x s) dWi, t > 0, i = 1, . . . , m, 
Jo j=l Jo 

where x G M. m and B\ A) G C 6 °°(]R m ) for alH = 1, . . . ,m and j = 1, . . . , d. We 
assume that the diffusion matrix A satisfies the following ellipticity condition 
at starting point x: 

(E) A(x)A(x)* > c*I m , for some c* > 0, where •* stays for matrix trans- 
position. Without loss of generality, we will suppose c* < 1. 

We recall that the first- variation process of X is the matrix- valued process 

OX^ 

^ i,j = ~dx ± ' i,j = l,...,m, 
which satisfies the following equation, written in matrix form: 

rt d rt 

dY t = I m + dB(X s )Y s ds + J2 9Ai{X s )Y s dW l s , 
Jo l=l Jo 

where dB and dAi are, respectively, the m x m matrices of components 
(dB)ij = djB 1 and {dA{)i j = djA\. By means of Ito's formula, one shows 
that Yt is invertible and that the inverse Zt := Y t _1 satisfies the equation 



(2.10) 



- j*Z 8 ^8B{X a ) - Y^dMXs)) 2 ^ ds 

V f z s dMx s )dw l s . 
1=1 J ° 



Additional notation. For k > 0, we define 



\ B \k = i + J2 Yl su p \ d « B \ x ) 



i=l 0<\a\<k X& 

(2.11) 



where \a\ is the length of the multi-index a. Then, for p > 1 and t > we set 
(2.12) e p (t):=e tP/2 ( tl/2 W^ A ^ P 
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and 

(2.13) e Z {t):=e t^H\B\i + \A\l) + \A\ l)r 

The constants in (2.12) and (2.13) naturally arise when estimating the mo- 
ments of the random variables Xt, Yt and Z%. Indeed, the results given in the 
following proposition can be easily obtained from (2.9) and (2.10) applying 
Burkholder's inequality and Gronwall's lemma. 



Proposition 2.2. For every p > 1 there exists a positive constant C v<m 
depending on p and m but not on the bounds on B and A and their deriva- 
tives such that, for every < s <t <T, 



(2.14) (i) E 

(2.15) (ii) 
for alli,j = l,. 



sup \X l r -X l s \ p 

L s<r<t 



<c p , m (t- s y/ 2 ((t-s)y 2 \B\o+\A\ y 



S upE[\(Z s ) lyj \P]<C p>m e^(t) 

s<t 



C p 



,171. 



For any t > 0, the iterated Malliavin derivative of Xt is the solution of 
a linear SDE. The coefficients of this equation are bounded, and hence it is 
once again a straightforward application of Gronwall's lemma to show that 
the random variables Dr^'^r^ Xt nave moments of any order which are finite 
and uniformly bounded in r\, . . . , r&. This is indeed the content of [16], The- 
orem 2.2.2. The following lemma highlights the explicit constants appearing 
in the estimates of the L p -norms of the iterated derivative, expressing them 
in terms of the bounds (2.11) on A and B. 



Lemma 2.1. For every k>l and every p > 1 there exist a positive inte- 
#er 7^ p and a positive constant Cfc jP depending on k,p but not on the bounds 
on B and A and their derivatives such that, for any t > 0, 

sup E[\Dil : ;;f k Xi\v} 

ri,...,r k <t 

(2.16) 

< CkJA^t 1 / 2 ^ + \A\ k )( k +^Pe p (tr^, 
for alii = l,...,m and (ji, ...,j k )e {1,. . . ,d} k . 

The proof of this result is based on some standard but rather cumbersome 
computations; hence we leave it for Appendix A.l. We rather give hereafter 
the proof of some estimates which follow easily from Lemma 2.1 and will be 
useful in the following sections. 
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Corollary 1. For any k > 1 and p > 1, i/iere exists a positive con- 
stant C k)P depending only on k and p such that, for any t>0, 

\2 



(2.17) 



(2.18) 



(i) E[\\D^Xl\f H ^ p < C Kp t k ' 2 \A\ k k -i{t ll2 \B\ k + \A\ k )^ 2 

xe p (ty*>v; 

(ii) ||^|| fcjP <C^<^(l + (tV^ /2 ) 



x\A\ k k _i(t 1/2 \B\ k + \A\ k )^ e p (t) k ^, 
where (i) holds for i = 1, . . . , m and (ii) for any (ft G C 00 (W ri ) . 

and Lemma 2.1, a sim- 



Proof. (i) Employing the definition of 
pie computation holds. 



n\D^xi\\iA x /v 



< Ck,p{ t 



^(p/2-1) 



E 



[0,*]* 



sup \D%;;;;gxlr 



\a\=k 



dr\--- dr k 



l/p 



< C k , p t k / 2 \A\ k k _i(t 1/2 \B\ k + \A\ k )^ 2 e p (t)^/P } 

\ k)P and write 

/ ft 
U(Xt)h, P = U[\<j>{X t )\*] + X}E[||2?W0(Jf, 



hence we get bound (2.17). 

(ii) We start from the definition of 



(2.19) 



<U\\ + ^E[\\D^ h U(X, 



h=l 



t/W^^h J 



Using the notation introduced in the proof Lemma 2.1, we have 

V 

D( h U(X t ) = Yl d kl ---d k J(X t )l[D^ d ^X k ', 

h,...,i v ={i,...,h} i=i 

where, with a slight abuse of notation, we have now written L>W for the 
generic derivative of order h. Repeatedly applying Holder's inequality for 
Sobolev norms and using bound (2.17), we get 

E[\\D( h U(X t )\\ P mh ] 



< c h , p Yl 

hi,...,hj/=l,...,h 

hi-\ \-h v =h 



E 
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<c KP u\\i E -p T[n\D^xi\^ hl ]V* 

i , , , i=l,...,m, -, rt 
hi,...,hij=l,...,h 1=1 

hiA \-h v =h 

By means of this bound, from (2.19) we get the desired estimate when setting 
Cfc,p > m a x {ch, p :h<k}. □ 

We need a last preliminary result on the inverse moments of the determi- 
nant of the Malliavin covariance matrix of Xt ■ This result is again achieved 
with some standard arguments, but, as in Lemma 2.1, the next lemma finds 
out the explicit constants appearing in the estimate of the L p -norms of 
det(o"x t ) _1 - 

Lemma 2.2. For every p> 1 and t > 0, 

(2.20) E[|detcjx t r p ] 1/p < C Pjm4 ef {mp+1) (tf^K m (t,c*), 
where 

K m (t,c*) = 1 + (± + ij + ^^{t^BUAf, + \\A\\\f^ 

for some positive constant C pmc i depending on p, m and d but not on the 
bounds on B and A and their derivatives. 

The proof is once again postponed to Appendix A.l. 

We now come to the main result of this section. We give an estimate of 
the L 2 -norm of the random variables H a involved in the integration by parts 
formula (2.8), when F = Xt- The proof follows the arguments of [6], proof of 
Lemma 4.11, but is given in the general setting of an integration by parts of 
order k G N, and moreover it takes advantage of the explicit bounds which 
have been obtained in Corollary 1 and Lemma 2.2. 

We give this result employing some slightly more compact notation, defining 

P k (t) = t 1 / 2 \B\ k + \A\ k , 
P A {t) = \A\ k P k+1 {t). 

Theorem 2.3. For every k > 1 there exists a positive constant C k = 
C k: m,d such that, for any multi-index a € {1, ... , m} k , any G G D°° and t> 0, 

\\H a (X t ,G)\\o : 2 

(2.21) <C k \\G\\ K2k+ i{r k ' 2 V t k ^ 1)/2 )(t m K m (t,c,)) k{k+3)/2 

x (P^(e8(t) V e 2 H 2 (i)) Ct (ef 2m+4 (t) V ef fc+4m+4 (t)) Cfc , 
where K m (t,c*) has been defined in Lemma 2.2, and 

cf) k = 3m(k + A) 2 . 
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Remark 2.1. Estimate (2.21) is rather involved. For our purposes, the 
most important elements are the dependence with respect to time of the 
factor i~ fc / 2 V t k i k ~ 1 )/ 2 anc [ the coefficient PS containing the bounds on 
the derivatives of the coefficients. We remark that the factor t m K m (t,c*) is 
bounded for t close to zero. Moreover, when t < 1, the factor t~ k l 2 \/t k ( k ~ 1 ^ 2 
reduces to t~ k / 2 . 

Proof of Theorem 2.3. We write a t = crx t f° r simplicity of nota- 
tion. We first use the continuity of 5 (see [6], Proposition 4.5) and Holder's 
inequalities for Sobolev norms to obtain 



(2.22) 



l#(ai,...,a fc )(^*'^9llo,2 

= \\ H {a k )( X t,H^ ait _ tak _^(X t ,G))\\ Qt2 
m 

J'=l 



0.2 



< Cm ||-g(a 1 ,...,oi fe _ 1 ) (X t , <?) || i,4 ^ || fa 1 )a k ,jhA DX t lll,8,«- 

J'=l 

To estimate the last factor we can directly use the definition of || • an d 
apply Corollary 1. The major part of the efforts in the rest of the proof will 
be targeted on the estimation of ||( cj r 1 )jj llfc,p- 

We claim that for any k > 1, p > 1 and for all i,j = l,...,m, 

(2 23) lite" 1 )*-*"** " Cfc ' p(t_1 Vt k ' 2 - l )(t m Km{t^)) l+k 

x P fc A (t)^ +2(fc+4)2 e P (t) Cfe -ef (mp+1) (t) c ^-, 

where 

^ = 2(fc + l)(m-l), 

and Cfc jP is a positive constant depending also on m, d but not on t and on 
the bounds on B and A and their derivatives. Iterating process (2.22) and 
repeatedly using estimates (2.23) and (2.17), one easily obtains the desired 
estimate 

\\H( ai! .., t a k )(Xt,G)\\ 0i2 

< Cfc ro rf||G|| fe 2 k+i(t m K m (t, c*)) fc 



J (t' 1 V t^ 2 - 1 )^ V t h ) l l 2 {t m K m {t, c*)) h 



k 
h=l 

x p A (i)^' fe +2(fc+4) 2 +(fc+i) 2 ) 

k 

x n e 2h+ 2(t) c ^ef (2h+2m+1) (t) c ^ 
h=l 
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<C k , m JG\\ k>2k+1 (t m K m (t,c^ k+3 y 2 

X (t- fc /2 vt ^l)/2 )p A (t) fe 

x (e 8 (t) V e 2k+ 2(t)f k ' m - d (e! 2m+ 4(t) V e^ +im+4 (t)f k ^ d . 

Proof of (2.23). We follow [6], proof of Lemma 4.11. We start from 
the definition of || • \\k >p and write 

/ k \ 1/p 

(2.24) \\(at 1 )i :j \\ k , p = fElK^-^^n + ^EtpW^- 1 )^-!!^]] . 

For the first term, we simply use Cramer's formula for matrix inversion, 

|K- 1 ) M | = (deta t )-Vp ) , 
j) minor of 

and bounds (2.17) and (2.20) and get 



where denotes the minor of at- We then apply Holder's inequality 



< ci 1 i l {E[det( ( Ti)- 2p ]EUp \\DX l t \\^ m ~ 1)p 



1 1/2 



(2-25) <c^^-P(t m ^ m ,(t,c,)r 

x e p (t) p -"^e 4(mp+1) (t) 

where K m (t, c*) is the constant defined in Lemma 2.2. To estimate the second 
term, as done in [6], proof of Lemma 4.11, we iterate the chain rule for D 

m 
a,b=l 

We take advantage of the notation introduced in the proof of Lemma 2.1 
and for (/?i, . . . , f3 k ) € {1, ... , m} k , k > 1, we write 

m 

hU-\Jl u ={l,—,k} ai '-> a " =1 
6i,...,6„=l 

(2.26) _ _ 

X (<T t )6„_ 1 ,a !/ (cr i )6„,j| 

x\D^(a t ) aiM ..-D^(a t ) av ,a 
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We repeatedly apply Holder's inequality for Sobolev norms to (2.26) and get 



<c 



u P 1 



(2) y- 

k,p,m / j 

- a a a . ft/JJ. X m m i i ■ rC 

feiH \-k v =k 



snp EHIK 1 )^ 1 
a,oi,...,ai/=l,...,m 
6,6i,.. .,&!/=!,... ,m 



(2.27) 



<c 



(2) 

k,p,m 



x sup ^ sup EO^- 1 )^!^ 1 ^ 
fei,...,fc !/ =l,... ; fe I a,6=l,...,m 

fclH h£v=fc 

v 

xj sup E[p( fe ')(a t ) aij6i | 

; =1 a ; ,6; = l,...,m. 



2 ( p nl/2' 



where, as in the proof of Corollary 1, we have written for the generic 

derivative of order fc/. To estimate D^ kl \at) ai ,bi we use bound (2.17) and get 

E[\\D^(a t )i 



\v 



< E 



( 2 '28) <=SE E [H D< " ,DX <*lll»«] 



£ 

/i=0 
A' 



D^D s Xl ■ D (k - h) D s x{ 



1/2 



h=0 



E[\\D( k - h ^DXf\\ 2 " 



h+li 



,1/2 



< cgt (fe/2+1)p i4* +2)p (t 1/2 i5| fc+1 + i^| fc+1 ) 2 ( fe+2 ) 2 %(t)^^-, 

where we have once again applied Holder's inequality for Sobolev norms in 
the second step. 

Using (2.28) together with (2.27), bound (2.25) and (2.24) and observing 
that t m K m {t) is greater than one for all the values of t, we finally obtain 

) (t 1 / 2 \B\ k+1 + \A\ k+1 /^ k +V 2 



x\A\f +k(k+2) ' 
xe P (t) Cfc ' p ef (mp+1) (t) Cfc -, 



for a positive constant Ck jP depending also on m,d. Estimate (2.23) follows. 

□ 
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2.4. Proofs of Theorems 2.1 and 2.2. We now come to the proof of the 
results stated in Section 2.1. We recall that an M m -valued random vector X 
is said to admit a density on an open set A G W 71 if Cx\a possesses a density, 
Cx being the law of X. It is equivalent to say that 

(2.29) nf(X)}= [ f{x)p{x)dx 

holds for all / G C&(R) such that supp(/) C A, for some positive p G L 1 (^4). 

We refer to the setting of Section 2.1. We recall that X = (X t ;t G [0,T]) 
denotes a strong solution of 

X* = x < + [ t b i (X s )ds + Y / f a){X s )dWl, 
Jo j=1 Jo 

(2.30) 

t G [0,T],i = 1, . . . ,m, 

where 6 and <r satisfy the assumptions (H1)-(H3). For k > 1 and / G C fc (M m ), 
we denote 

(2.31) I/U,BkM = 1 + E SU P M x )l 

\a\<k xeB *to°) 

where the sum is taken over all the multi-index a G {1, . . . , m} k . Let us define 
the following "local" version of the constants appearing in the estimates of 
the previous section: 

Pk{t,Vo) = tl/2 \ b \k,B 5R (y ) + W\k,B SR (y ), 

P k(t,yo) = \a\ kjBsM P k+ i(t,y ), 

p i(t,yo)=t 1/2 (\b\ ltB5R{yo) + \o-\1 BbrM ) + \a | l,B BR (y ), 

|2 ^ 

C m (t,y ) = t m + -5(^(1 + P%(t,y )), 
c yo 

e p (t,y ) = eMi p/2 Pi(t,yo) p ), 

e z p (t,y ) = eMt p/2 Pi Z (t,yo) p )- 
In order to prove Theorem 2.1, we simplify this rather heavy notation in- 
troducing a constant that contains the factors appearing in estimate (2.21) 
in Theorem 2.3 (recall the constant 4>k defined there) 

e fc (i,y ,7) = C m (t,yo) mfc(mfc+3)/2 ^ fc (t,yo)^ fc+(mfc+2)2 
x (e&(t,yo) Ve 2 mfc+ 2 (t,yo)) 7 
x (ef2 m +4(*>yo) Vef mfe+4m+4 (t,y )) 7 - 



D C/, \ /j.1/2 1 7i I i3 I 1 2 \2(m+l) 



16 S. DE MARCO 

As addressed in Section 2.1, the following theorem is a more detailed 
version of Theorem 2.1. In particular, it provides the explicit expression of 
the constant A k appearing in estimates (2.2) and (2.3). 

Theorem 2.4. Assume (HI), (H2) and (H3). Then, for any initial con- 
dition x £ R m and any < t < T , the random vector X± admits an infinitely 
differentiable density pt, yo on Bji(yo). Furthermore, for every k>l there 
exists a positive constant C k = C kiiriiC [ such that, setting 

(2.32) A k (t, y ) = C k R- mk (P (t, y ) mk + @ k (t, y , C k )) 
and 

P t (y) = P(inf{|A s - y\ :s £ [(* - 1) V t/2,t}} < 3R), 
then one has 

(2.33) p tm (y) < P t (y ) (l + ^ J A 3 (t A 1, y ) 

for every y £ Br{i)q). Analogously, for any a E {1, . . . ,m} fc , > 1, 

(2.34) \d a Pt,y (y)\ < Pt(yo)(l + — ^) A 2fc+ 3(t A l,y ) 

/or every y£B R (y ). 

To prove this result we rely on the following classical criterion for smooth- 
ness of laws based on a Fourier transform argument (cf. [16], Lemma 2.1.5). 



Proposition 2.3. Let fi be a probability law on R m , and /!(£) = f R e i{ £>v) x 
n(dy) its characteristic function. Iffi is integrable, then [i is absolutely con- 
tinuous w.r.t. the Lebesgue measure, and 

(2-35) p(y) = -i^ / e-^MCO* 

is a continuous version of its density. If moreover 

(2.36) / \d\ k \m\dZ<oo 

holds for any k £ N, then p is of class C°° and for any multi-index a = 
(ai,.. .,a k ) e{l,...,m} k , 

d aP (y) = (~i) k jf (n^ j ^' V) m 
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Proof of Theorem 2.4. Step 1 ("localized" characteristic function) . 
Fix a t in (0,T). Let <f> R G C^°(M m ) be such that l Bfl(0) < fl < l BaH(0) 
and |</>hU < 2 k R~ k . We first observe that if uiq = K[^>n(Xt — yo)] is zero, 
then it just follows that p = is a density for on B R (yo). Otherwise, we 
consider Ct,y the law on R m such that 

(2.37) / f(y)£ t , yo (dy) = —E[f(X t )<f> R .(X t - y )}, 

for all / G Cf,(]R m ). If £j >3/0 possesses a density, say Pt tVo , it follows that 
Pt, yo (y) ■= m oPt, yo is a density for X t on B R (y ). Indeed, for any / G C b such 
that supp(/) C B R (yo), (2.37) implies 



f{y)Pt, yo {y)dy = / f{y)m p^ yo {y)dy 



= m f(y)£ t , yo {dy) 

JR m 

=nf(Xt)]. 

If the characteristic function of Ct,y 

Pt,w(0 = / e^Ct^dy) = —Eie^MXt - yo)] 

is integrable, then by Proposition 2.3 Ct >yo admits a density. Hence, we focus 
on the integrability of Pt,y ', in particular, we show that condition (2.36) of 
Proposition 2.3 holds true for all k G N. 

Moreover, the inversion formula (2.35) yields the representation for pt yo 

Pt,yo (V) : = m 0Pt,y (y) = 77TT^ / e~^' V) Pt,y (£) <% 

(2.38) 

= TT^ I e- i{U W { *' Xt) MXt ~ yo)] dC- 
(2vr) m J Rm 

Step 2 (localization). We define the coefficients 

V(y) = b\iP(y-y )), 

(2.39) 

°*j(v) = <TjMv-Vo)), 
where ip G C oc (M. m ;M m ) (a truncation function) is defined by 

(y, if|y|<4fl, 

m = \ 5 y ) if \ y \>5R, 

y \y\ 

and ip(y) G B§ R (0) for all y G M. m . t/> can be defined in such a way that, 
for all i = l,...,m, ||^||i < 1 and ||^|L < 2 k ~ 2 R~ { - k ~ 1 ^ for all fe > 2. As 
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a consequence of (HI), the b and a denned in this way are C£°-extensions 
of b\ 

B 4R (y ) an d a \B 4R (y )- Furthermore, there exist constants c = c\-^ m such 

that 

\b Z \k<Ck,rnR \b l \k,B 5R (y ), 

(2.40) 

and by (H2), for any y G B 3R (yo) the matrix a(y) is elliptic 
(2.41) oo*{y) > c yo , R Im, V G B 3R {y ). 

For y G M m we denote by X(y) = (X s (y);0 < s <t) the unique strong solu- 
tion of the equation 

Xi(y) = y l + [ S b l (X u (y))du + Y] f a){X u {y))dWl 
Jo j=1 Jo 

(2.42) 

< s < t,i = 1, . . . ,771. 

Let now < S <t/2 A 1. We employ an up-down crossing argument to 
estimate the increments of X in the neighborhood of yo by replacing them 
with the increments of X. More precisely, let v = Vt8 an d t = Tts be the 
stopping times defined by 

u t s = inf{s > t - 5 : X s G B 3R (y )}, 

(2.43) 

T t) s = inf{s > i/ t|<5 : X s <£ B AR (y Q )} 

and inf{0} = oo. Suppose that 4> R (X t — yo) > 0, so that X t G B2 R (yo) and 
v < t. On this set, if v > t — 5, then |^ At — X v \ > ii. This implies 
\X t /\ T -y{Xy) — X v \ = \X t /\ T — X u \ > R. Here we are employing the fact that 
on the interval [v, r], X stays in B^ R (yQ), hence in the region where the trun- 
cated coefficients b, a coincide with the original ones b, a. On this interval, 
both X and X satisfy (2.42) for which pathwise uniqueness holds; hence we 
can replace X by X and employ the flow property for X. Notice that flow 
property may not hold true for X [due to possible lack of uniqueness for the 
couple (b,a)], but it always does for X. 

Analogously, if v = t — 5 and r<t, then \X T — X v \ = \X T - V (X U ) — X„\ >R. 
In both cases, sup 0<s<s \X S (X U ) — X„\ > R. Hence, we conclude that 

{<f> R (X t - yo) > 0} = {4> R {X t -y )>0,t-5 = u<t<T} 

u \Mx t - yo) > o, sup \x s (x u ) - x v \ > r\ 

1 0<s<5 J 

and pt,y rewrites as 

m pt, yo (0 = E[e l ^MX t - yo)l {fe(Xi „ w)>0 ,sup < s < 5 \x a {x v )-x v \> R }\ 
+ E[e i{ t' Xt ^ R (X t - yo)l{4, R {X t -y )>0,t-8=u<t<T}}- 



SMOOTHNESS OF DENSITIES AND ASYMPTOTIC ESTIMATES 



19 



We now claim that for all q > the following estimate holds: 



<p R {X t - y ) > 0, sup \X S (X U ) -X U \>R 

v 0<s<5 

(2.44) 

< C( ,, mJ R-^/ 2 Po(5,yo)"IP( pf \X s -y \<3R 

\t—S<s<t 

for some positive constant c q ^ m . Estimate (2.44) will be proved later on. 
On the other hand, 

W^ Xt H R {X t - Ito)l WH (Xt- OT )>0,t-«=l/<t<T}] 

(2.45) = ElEie^W <t> R (Xs(y) - yo)\X t - S = y]l {t ^ =u<t<T} ] 

< n\X t -8 - Vol < 3-R) sup \E[e^M)^ R (X 5 (y) - y Q )}\. 
yeB 3R (y ) 

Step 3 (integration by parts). We apply integration by parts formula (2.8) 
to estimate the last term in (2.45). By (2.40), (2.41) and Lemma 2.2, X$(y) 
is a smooth and nondegenerate random vector for any 5 > and y € B 3R (yQ). 
Then, for a given k > 1 we define the multi-index 

a = (!,...,!,..., m,...,m) , 

k times k times 

such that \a\ = km. Hence, recalling that d Xk e l ^' x ^ = i^ k e l ^' x \ 

|E[e i «^>^(X 5 (y)-y )]| 
(2-46) < 1 \E[d a e^^MX 5 (y)-y )]\ 

Ui=i\C\ 

< Y^j^E[\H a (X s (y), ( f )R (X s (y)-y ))\}, 

for any y G B 3R (y ). _ 

We need to separately estimate \\<f>(Xs(y) — yo)\\\ a \ 2H+ 1 - By Corollary 1, 
this is given by 

U(Xs(y) - yo)\\ mk , 2 ^ < c k , m R- mk (i + 5 1/2 )k|™t lii?5flfeo) 

xP m k(yo,S) (mk+2)2 e 2mk+ i(5y^ 

< c k , m R- mk P^ k (y , «5)( mfc+2 ) 2 e 2mfc+1 (5) c <- 

for some positive constant c ktTn . Then, from (2.44), (2.45), (2.46) and The- 
orem 2.3 it follows that 

(2.47) mo|Pt,yo(£)l < C k , q P R (6,t,y )I k , q (t,5,y ) 
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for some constant C ktQ depending also on m and d, with 
P R (5,t,y )=F( inf \X s -y \<3R) 

\t—S<s<t / 

and 

Ik, q (M,Vo) = R- q S^ 2 P (S,y y + „ ™ k 5- mk / 2 e k (5,y ,C k , q ). 

lli=i Is I 

Estimate (2.47) holds simultaneously for any £ G W' 1 , < S < t/2 A 1, g > 
and A; > 1. The constant 8fc(<5, yo 5 Cfc,<j) appears when applying estimate (2.21) 

Step 4 {optimization). We show that for any £ and any I > 1, <5 can always 
be chosen in such a way that there exist q and fc such that I kj q(£,5,yo) goes 
to zero at oo faster than (n^il?i)" (/+2) • 

Denoting = II™ ilf I, we set 

( 5:=,5(0=t/2AlA||eir a 
for some a > that is to be identified hereafter. For this choice of 5, 

PR(6(Q,t,yo) <f( inf \X S - y \ < 3r) = P t (y ) 

\t/2V(t-l)<s<t / 

and 

h, q & m,Vo) < R- q (U\r qa/2 A (t A l) g / 2 )P (t A 1, y ) 9 

(2.48) + R-mk^-k(l-ma/2) y ||£||-k(t A !)-mfc/2) 

x & k (tAl,y ,C kt q), 

since 5 — > Pq(8, yo) and 5 — > & k (S, yo,C k)Q ) are increasing; hence Po(S(t;),yo) < 
Po(t A l,yo) and the same holds for Q k . 

We consider the leading terms determining the decay of I k ,q(£, 2/0 ) 
with respect to £ and impose 

qa , / ma\ 

(2.49) \=\ i --y 

Setting a = 1/m, (2.49) yields g = mfc, hence f = jfe(l - *f) = §. Therefore, 
we get the bound 

Ik, qi (ZM),Vo) < R- mk (M\\- k/2 A (t A l) mk / 2 )P()(t A 1, yo) mk 

(2.50) + ir mfe (H£ir fc / 2 v (t a i)- mfc/2 iieir fc ) 

x Q k (t A l,yo,C fci9 *) 

with gjf; = mA;. Estimate (2.50) holds for any k > 1 and £ 7^ 0, and then it 
proves that the function pt,y (y) defined in (2.38) is in fact well defined and 
infinitely differentiable with respect to y. 
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Let us come to estimate (2.33). We take (2.38) and cut off the integration 
over a region / of finite Lebesgue measure on which ||£|| = nULil^l remains 
smaller than a given constant. That is, we write 

1 



Pt,y {y) 



< 



(2tt)' 
1 



< 



(2tt)' 



E[Mx t -yo)]dt; 

5 -*«.f>E[ e *«.*>0 fl (x t -ito)]de 

P(|X t -y |<2 J R)A m (/) 



+ 



where A m denotes the Lebesgue measure on M m . As we have seen, the last 
term is such that 



i' 



<R- mk P (tAl,y ) mk [ |er fc/2 ^ 

+ R- mk Q k {t A 1, y , f (< A l)~ mfc / 2 jf 



+ 



c n{€:|€|<(*Ai)-»} 

lf|- fc/3 «Y 

/ c n{^: |?|>(tAl)- m } / 



Now, since 



/ ir fc/2 #< / ir* /a de=4 1) <°° 

Jl c r\{£, ■. \£\>(tf\i)- m } Jl c 



and 



(tAl) 



-mfc/2 



|Cr fc ^<(*Al)- mfc / 2 4 2) <oo 



'/ c n{M£|<(tAi)- m } 
hold for any k > 3, we then take = 3 and get the estimate 

Pt m (V) < C*P t (y )[l + R- 3m (Po(tAl,yof m + (t A l)- 3m / 2 G 3 (i A 1, y , C m4 ))} 

for every y G Bn(yo), for a positive constant C*, estimate (2.33) then follows. 
For estimate (2.34) on the derivatives we proceed in the same way, observing 

that for a € {1, ... ,m} 1 , |£|~ fc ^ 2 x Ylj=i \C° J \ 1S integrable at oo as soon as 
k > 21 + 3. 



22 S. DE MARCO 

Proof of (2.44). We remark that {(p R (X t - y )} ^{t-5<v<t}C 
{t-5<v<t,X u £ B 3R (y )}, hence 

<> R (X T - yo) > 0, sup \X S {X U ) -X V \>R 

0<s<<5 

<p(t - 5 <v <t,X u eB 3R (y ), sup \X S {X U )-X U \>R 

v 0<s<<5 

< R~ q F(t -8<v<t) _sup E [ sup \X s {y) - y\ q 

Using boundedness of coefficients of (2.42), it is easy to show that 
E[ sup \X s {y)-y\i] < c qjm 5^ 2 P (5, y Q f 

l 0<s<8 J 

for some positive constant c g>m . □ 

As addressed in Section 2.1, the constants appearing in the definition of 
Afc (2.32) can be considerably simplified under assumptions (Hi') and (H4), 
resulting in some polynomial-type bounds. The following result corresponds 
to Theorem 2.2; in the presents statement, we explicitly give the expression 
of the exponent q' k {q) appearing in bound (2.6). 

Theorem 2.5. Assume (HI') and (H3). 

(a) For any initial condition x £ M. m and any < t < T , Xf admits a smooth 
density on M. m \B v (0). 

(b) Assume (H4) as well. Then the constant defined in Theorem 2.4 is 
such that 

(2.51) A k (t,y )<C ktT (l + \y \^ q) ), 

for every < t < T and every \y\ > rj + 5. The exponent q' k {q) is worth 
q' k (q)=mk(q + 4)(m+l)(mk + 3) + 2qm((j) mk + (mk + 2) 2 ). 

(c) If moreover sup < s <t \X S \ has finite moments of all orders, then for every 
p > and every k>\ there exist positive constants Ck,p,T such that 

b^)l<^3,p,T(l + ^W, 

(2-52) i 

\d a p t (y)\ < C k , p , T (l + tm{2k+3)/2 j \y\~ p , a € {1, ... , m} k , 

for every < t < T and every \y\ > rj + 5. 

The C ktPt T are positive constants depending also on m,d and on bounds (2.4) 
and (2.5) on the coefficients. 

Proof, (a) We no longer need to distinguish between yo and the (close) 
point y where the density is evaluated; hence we just set y = yo and consider 
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suitable radii. For \y\ > 77, we set R y = jq dist(y, B v (0)) A 1. By (HI'), b and a 
are of class on B^R y {y) and satisfy (H2) on B^R y {y). From Theorem 2.4 
it follows that Xt admits a smooth density on Bji y (y). This holds true for 
every ball Bn y (y) with center y in M. m \ 5^(0); hence statement (a) follows. 

(b) Without loss of generality, we take R = 1. As a consequence of (2.4), 
the constants introduced before Theorem 2.4 can be bounded as follows, for 
0<t<T and \y\ > 77 + 5: 

Pk(t,y)<c£ ) (l + (\y\+5)i)<c£ ) \y\i, 

PZ{t,y)VP?{t,y)<cV\y\ 2q i 
(2.53) P%(t,y) < c«(| y | 4 + \y\ 2 ) 2(m+1) < c«|y| 8 ( m+1 \ 

^(^)<^TT)l2/l Wl) l^l 8(m+1) 

< c (1) |y| 2(5+4)(m+1) , 

for some constants and cj^ depending also on m, q and on the bounds 
on b,a and their derivatives in (2.4) and (2.5). 

The exponential factors e and e z must be treated on a specific basis. 
Indeed, e.(t,y) and e z (t,y) may explode when \y\ — > +00. Nevertheless, ex- 
plosion can be avoided stepping further into the optimization procedure set 
up in the proof of Theorem 2.4. More precisely, we restart from step 4 and 
force the state variable y to appear in the choice of 5, setting 

^,y)=t/2AlA|er a A|y|- 4 ". 

Now, whatever the value of p is, e p and e z are reduced to 

e p (5(Z,y))<e p (lA\y\- 4 i,y) 

< exp(4 1) (l A \y\~ 2qp )(l + 1 A \y\- 2q ) p \y\ qp ) 
<exp(c p (lA|yr 2 ^)|y|^) 

< exp (c p ) 

and 

e z p {5{t,y))<e z p {lf\\y\-\y) 

= exp (cS 1} (l A \y\~ 2qp )(l + 1 A \y\~ 2q ) p \y\ 2qp ) 

< exp(cp). 

We then perform the integration over £ as done in the last step of the proof 
of Theorem 2.4, and employing (2.53) we obtain estimate (2.51) for A*., 
for \y\ > T] + 5. The value of q' k (q) is obtained from the definition of 0& 
and (2.53). 
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(c) From boundedness of moments of sup s<t \X S \ , for any interval /( C [0, t] 
one can easily deduce the estimate 

P(inf{|X s -y|:sG/ t }<3) 

(2.54) <P(sup{|X a |:a€l t }>M-3) 



< 71 i 1 x E 

" (M-3) r 



sup \ X S \ 

s<t 



\y\ r 



for any r > 0, < i < T and \y\ > 3. It is then easy to obtain the desired 
estimate on pt with Theorem 2.4: for a given p > 0, we employ (2.33) with 
yo = U and (2.54) with r > p + q' 3 (q). Similarly, to obtain the estimate on 
derivatives, one employs (2.34) and (2.54) with r > p + q' 2k+s (q) . □ 

3. A square root-like (CIR/CEV) process with local coefficients. We ap- 
ply our results to the solution of (1.1). We will be able to refine the polyno- 
mial estimate on the density at +oo giving exponential-type upper bounds. 
Under some additional assumptions on the coefficients, we also study the 
asymptotic behavior of the density at zero, that is, the point where the 
diffusion coefficient is singular. 

We first collect some basic facts concerning the solution of (1.1). Let us 
recall the SDE 



(3.1) 



dX t = (a(X t ) - b(X t )X t ) dt + j(X t )X? dW t , t>0,ae [1/2, 1), 
X = x > 0. 



When a = 1/2 and a, b and 7 are constant, the solution to (3.1) is the cele- 
brated Cox-Ingersoll-Ross process (see [5]), appearing in finance as a model 
for short interest rates. It is well known that, in spite of the lack of globally 
Lipschitz-continuous coefficients, existence and uniqueness of strong solu- 
tions hold for the equation of a CIR process. If a > 0, the solution stays 
a.s. in = [0,+oo); furthermore, a solution starting at x > stays a.s. in 
M> = (0, +00) if the Feller condition 2a > 7 2 is achieved (cf. [14] for details). 
The following proposition gives the (straightforward) generalization of the 
previous statements to the case of coefficients a, b, 7 that are functions of 
the underlying process. The proof is left to Appendix A. 2. 



Proposition 3.1. Assume: 

(sO) a G [1/2,1); a, b and 7 G with a(0) > and -f(x) 2 > for every 
x>0. 

Then, for any initial condition x > there exists a unique strong solution 
to (3.1) which is such that ¥(X t > 0;t > 0) = 1. Let then x > and tq = 
infji > : X t = 0}, with inf{0} = 00. 
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• If a > 1/2 and 

(si)' a(0) > and z h-» ^^^ia-i integrable at + , 
i/ien 

(3.2) P(r =oo) = l. 

• If a = 1/2 and 

(si) -if is integrable at zero, 

(s2) i/iere exists x > suc/i £ftai > 1 /or < x < x, 
then the same conclusion on tq holds. 

When X is a CIR process, the moment-generating function of Xt can 
be computed explicitly, leading to the knowledge of the density. Setting 
Lf = (1 — e _b ')7 2 /46, then Xt/ 'Lt follows a noncentral chi-square law with 
5 = 4a/7 2 degrees of freedom and parameter Q = Axb / {e bt — 1)) (recall 
that x is here the initial condition). The density of Xt is then given by 
(cf. [14]) 

Pt (v) - ^e-y'^) (v\ 5/2-1 t {y/{4Lt)r c v > o 

PM>- 2 S/2 Lt e { L J Z, n!r(J/2 + n )^' y> °- 

We incidentally remark that pt is in general unbounded, since y^ 2-1 diverges 
at zero when 5/2 — 1 = 2a/ r y 2 — 1 is negative (in fact, fixed a value of 5/2 — 1, 
there exists a n > such that JprP* is unbounded) . 

The standard techniques of Malliavin calculus cannot be directly applied 
to study the existence of a smooth density for the solution of (3.1), as the 
diffusion coefficient in general is not (depending on 7) globally Lipschitz con- 
tinuous. Actually, Alos and Ewald [1] have shown that if X is CIR process, 
then Xt, t > 0, belongs to D 1 ' 2 when the Feller condition 2a > 7 2 is achieved. 
Higher order of differentiability (in the Malliavin sense) can be proven, re- 
quiring a stronger condition on a and 7, and the authors apply these results 
to option pricing within the Heston model. If we are interested in density 
estimation, the results of the previous sections allow us to overcome the 
problems related to the singular behavior of the diffusion coefficient and to 
directly establish the existence of a smooth density, independently from any 
Feller-type condition [provided that (sO) is satisfied]. More precisely, we can 
give the following preliminary result: 

Proposition 3.2 (Preliminary result). Assume (sO) and let 0,6,7 be 
of class C£° . Let X = {Xt,t > 0) be the strong solution of (3.1) starting at 
x > 0. For any t > 0, X t admits a smooth density p t on (0, +00). pt is such 
that lim^oo p t (y)y p = for any p>0. 

Proof. It is easy so see that, under the current assumptions, the drift 
and diffusions coefficients of (3.1) satisfy (HI') with i] = and (H4) with 
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q = 1- (H3) holds as well, by Proposition 3.1. As the coefficients have sub- 
linear growth, for any t > 0, sup s<t X s has finite moments of any order. The 
conclusion follows from Theorem 2.5(c). □ 

3.1. Exponential decay at oo. In order to further develop our study of 
the density, we could take advantage of some of the generalized-chaining 
tools settled by Viens and Vizcarra in [17]. In particular notice that, in 
order to estimate the density by means of Theorem 2.5, we need to deal 
with the probability term Pt(y) appearing therein. For our present purposes, 
we can rely on alternative strategies involving time-change arguments and 
the existence of quadratic exponential moments for suprema of Brownian 
motions (Fernique's theorem) in the current section, and a detailed analysis 
of negative moments of the process X in Section 3.2. 

From now on, condition (sO) is assumed, the coefficients a, b, 7 are of 
class C£° and (X t ;t < T) denotes the unique strong of (3.1) on [0,T], T > 0. 
We make explicit the dependence with respect to the initial condition de- 
noting pt(x, ■) the density at time t of X starting at x > 0. The following 
result improves Proposition 3.2 in the estimate of the density for y — > 00. 

Proposition 3.3. Assume that 



for every y>x + l, with C = 2 3 ~ 2a + 2|-y|§(l - a) 2 . The C k {T) also depend 
on a and on the coefficients a,b and 7. 

Remark 3.1. In the case of constant coefficients and a = b = 0, the 
bound (3.4) can be compared to the density of the CEV process as pro- 
vided, for example, in [7], Theorem 1.6 (see also the references therein). The 
comparison shows that our estimate is in the good range on the log-scale. 

Proof of Proposition 3.3. In the spirit of Lamperti's change-of-scale 
argument (cf. [12], page 294), let tp G C 2 ((0,oo)) be defined by 



(3.3) lim b(x)x 1 - a > -00. 



Then there exist positive constants 70 and C^iT), k>3, such that 




and 





h\ov 



1 



,tt 



dy 



Mo(l -a) 



1 



x 



l-a 
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so that ip'(x) = | 7 | ox c ■ Let moreover 6 £ C£°(M) be such that l[2,oo)^#^l[i,oo) 
and 9'<l.We set 



p(x) 



0(x)ip(x), x > 0, 
0, x<0, 



so that p is of class C 2 (M). We define the auxiliary process Yt = Xt — x, 
t > 0, which is such that P(lt > —x) = 1, t > 0. An application of Ito's 
formula yields 

ft 



where 



and 



p{Y t ) = [ f(Y s )ds + M t 
Jo 

f(y) = p'(y)(a(y) - b(y)y) + y'(yh(y) 2 y 2a 

M t = f p'(Y s )j(Y s )Y?dW x . 
Jo 



The key point is the fact that / is bounded from above on (— x,oo) and, 
on the other hand, M is a martingale with bounded quadratic variation. 
Indeed, / is continuous, it is zero for y < 1 and for y > 2 one has 

I[y) h\oy a h\o y 2 | 7 | y ' 

and hence, recalling that a is bounded, lim 2/ _ s . 0O f(y) < oo is ensured by 
condition (3.3). Then we set C\ = sup y>0 f(y). For M, one has 

(M) t = I V(y,)W^A« 

Jo 

(3.6) < A^ + ^^y,))^^) 2 ^ 20 ^ 

Jo 

<2(<^(2) 2 + l)t, 

and hence we set C2 = 2((p(2) 2 + 1) = 2( . 2 a!!j*\ 2 + !)• Now, since p is strictly 

I T 1 V a ) 

increasing, {Yt > y} = {p(Yt) > p(y)} for any y > 0. Moreover, {p(Y t ) > 
p(y)} C {M t + Cii > p(y)} C {2M 2 + 2C 2 t 2 > p(y) 2 }. We set I t = [(t - 
1) V t/2,t] and r = inf{s > 0:Y S > 3/2}. The quadratic variation of M is 
strictly increasing after r, since (p(Y T \/ t )^(Y T \/t)Y°^ t ) 2 > 0: hence, by Du- 
bins and Schwarz's theorem (cf. Theorem 3.4.6 in [12]) there exists a one- 
dimensional Brownian motion (bt;t > 0) such that M rV t = &(M) TV r Clearly 
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one has {Yt > 2} C {t < t}, so that for y > 2 

P t (y) = ¥(3s €l t :Y s >y)< F(3s G I t :Y s > y,r < s) 

< F(3s e I t :2M 2 + 2Cfs 2 > p{y) 2 ,r< s) 

< F( sup (2M 2 + 2C 1 V) > p(yf 

V<s<* 



<P( sup 6 2 A/)s +C 2 t 2 >ip(y) 2 

1 



<P sup b 2 + C 2 t 2 >-p(y) 2 

\s<C 2 t z 

We now employ the scaling property for the Brownian motion (b s ; s > 0) ~ 
(yfab s / a ;s > 0), a > 0, and Fernique's Theorem (cf. [11], page 402). The latter 
tells that there exists a positive constant 70 such that exp(7o sup s<1 b 2 ) is 
integrable, hence 

(C 2 
l0 supb 2 s + l0 -±t> 7 ^-p(y)' 
s<l ^2 

< exp ( - 7 o4^1 E[e^ 2 / C2 ^ su P*<i b *] 

< Coexp (- 7o ^! +7 „| t 

where Co = E[exp( 7 o sup s<1 6 2 )] is a universal constant. The estimates on 
the density of Xt and its derivatives now follow from Theorem 2.4 [esti- 
mates (2.33) and (2.34)] and Theorem 2.5(b), using Xt — x = Yt, the value 
of the constant C2 and taking, for example, R = 1/6. □ 

3.2. Asymptotics at 0. We have established conditions under which the 
solution of (3.1) admits a smooth density pt on (0, +00). According to Propo- 
sition 3.1, the process remains almost surely in M + : this trivially means that 
for any t > 0, Xt has an identically zero density on (— 00, 0), which can be 
extended to when tq = 00 a.s. We are now wondering what are sufficient 
conditions for pt to converge to zero at the origin, hence providing the exis- 
tence of a continuous (eventually differ entiable, eventually C°°) density on 
the whole real line. 

What we have in mind is the application of Theorem 2.5 to the inversed 
process Yt = (considered on the event {tq = 00}). An application of Ito's 
formula yields 

(3.7) dY t = J a (Y t ) dt - Wt)Y t 2 - a dW t , 
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where 

J a {Y t ) = -a{Y t )Y t 2 + b(Y t )Y t + j(Y t ) 2 Y?- 2a 

with the notation f(y) = f(l/y), y > 0, for / = a,b,a. Equation (3.7) has 
super-linear coefficients, in particular condition (2.4) of Theorem 2.2 holds 
with q = 2. Willing to apply Theorem 2.5(c), we first need some preliminary 
results on the moments of Y. The proof of the next statement is based on 
the techniques employed in [3], proof of Lemma 2.1, that we adapt to our 
framework. 



Lemma 3.1. 



(1) If a > 1/2, assume (si)'. Then for any initial condition x>0, for any 
t>0 and p>0 

1 " 



(3.8) 



E 



SU PyP 

s<t As 



<c, 



for some positive constant C depending on x,p,a,t and on the coeffi- 
cients of (3.1). 
(2) If a = 1/2, then assume (si) and (s2) and let 



(3.9) 



,* ,. 2a(a?) , 
1* = lim ^-nr > 1. 



X-*0 



7(x) 2 



Then (3.8) holds for p > such that 
(3.10) p+Kl*. 

PROOF. Let r n be the stopping time defined by r n = inf{t > 0: X t < 
1/n}. The application of Ito's formula to Xf Ar , p > 0, yields 

ftf\T n 



(3.11) 
where 



1 









E 



b(x) p (p+l, x 2 „2o-l 



=p 1 -r 

Yy ' y X P XP+ 1 



■y(x) x — a{x) 



x > 0. 



It is easy to see that, if 
(3.12) 



lim g(x) < 0, 



there exists a positive constant C such that (p(x) < p^-y + C, for every a; > 0. 
If (3.12) holds, from (3.11) we get 



E 



1 

sup^p — 

s<t A sAt„ 



<^ + Ct + p\b\oj Q E 





1 


f E 


u<s a uAt„. 


Jo 
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and hence, by Gronwall's lemma, 

1 



(3.13) 



E 



sup 

a<t 



X?, 



< 



— + ct 

xP 



J>\b\ot 



We verify (3.12), distinguishing the two cases. 

Case a > 1/2. We simply observe that \im x ^og(x) = — a(0) < 0. Esti- 
mate (3.8) then follows by taking the limit n — > oo in (3.13) and using Propo- 
sition 3.1 under assumption (si)'. 

Case a = 1/2. We have g(x) = P ^-^j 2 {x) — a(x). If p satisfies (3.10), 
then (3.9) ensures that \im x ^Qg(x) < 0. We conclude again taking the limit 
n — > oo and using Proposition 3.1 under assumptions (si) and (s2). □ 

We are now provided with the tools to prove the following: 

Proposition 3.4. 

(1) If a > 1/2, assume (si)'. Then for every t > 0, every p>0 and every 
k > the density pt of Xf on (0, +oo) is such that 



(3.14) 



lim . y p \p t (y)\ =0, 
lim .y- p \pl k \y)\ = <>- 



(2) If a = 1/2, then assume (si) and (s2) and define I* as in Lemma 3.1. If 
(3.15) T> 3 + ^(2) 

[where q'(-) has been defined in Theorem 2.5], then 
(3.16) 



for every <p <l* 
(3.17) 

then 
(3.18) 

for every 0<p<l* 



lim y p p t {y) = 

(3 + ^(2)). Moreover, if 
l*>2k + S + q' 2k+3 (2), 

lim -y- p \ P ?Hy)\=0 
(2k + 3 + q' 2k+3 (2)). 



Proof. We apply Theorem 2.2 to Y = 1/X . For simplicity of notation, 
we write p for pt and py for the density of Yf. As Y satisfies (3.7), from 
Theorem 2.2(b) it follows that the bound (2.51) on holds with q' k {q) = 
q' k (2). Hence, from Theorem 2.2(c) it follows that 



(3.19) 



lim p Y (y>)\y'\P' =0 
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and, for a given k > 0, 

(3.20) lim pf\y')\y'\P' =0, 

if sup s<t Y s has a finite moment of order r > p' + q^k+s^)- 
Now, it is easy to see that 

(3.21) p(y) = \pY[-), 

y \yj 

and hence, after some rather straightforward computations, 

/ 1 \ 2(fc+i) k i jv / 1 \ 

(3.22) |^%)|<C fc Qj EE^Q)< 0<»<1- 

Once again, we distinguish the two cases. 

Case a > 1/2. If 1/2 < a < 1, by Lemma 3.1, (3.19) and (3.20) hold for 
any p' > 0. Then (3.14) easily follows from (3.21) and (3.22). 

Case a = 1/2. By Lemma 3.1, (3.15) is the condition for sup s<t Y s to have 
finite moment of order strictly greater than 2 + q' 3 (2) + p, with p < I* — (3 + 
q' 3 (2)). By (3.21), in this case (3.19) holds true with k = and p' = 2 +p, 
and hence (3.16) holds. Similarly, by (3.22) estimate (3.18) holds if (3.19) 
holds with p' = p + 2(fe + 1). The latter condition is achieved if sup s<t Y s has 
finite moment of order strictly greater than 2(k + 1) + 9 2 (fc+3)(^) ~^~Pi whidi 
is in turn ensured by (3.17). □ 

Remark 3.2. Proposition 3.3 states that pt decays exponentially at 
infinity for any value of a, as far as condition (3.3) holds true. When a > 1/2, 
Proposition 3.4 states that pt and all its derivatives tend to zero at the origin, 
while the price to pay for the same conclusion to hold is higher when a = 1/2 
[cf. conditions (3.15) and (3.17), which become rapidly strong for growing 
values of k\. With regard to this behavior at zero, we recall that Proposi- 
tion 3.4 only provides sufficient conditions for estimates (3.16) and (3.18) to 
hold. We do not give any conclusion on the behavior of the density at zero 
when condition (3.15) [or (3.17) for the derivatives] fail to hold. 

APPENDIX 

We collect here the proofs of some of the more technical results. 
A.l. Proofs of Lemmas 2.1 and 2.2. 

Proof of Lemma 2.1. We refer to the notation introduced in the proof 
of [16], Theorem 2.2.2, allowing us to write the equation satisfied by the 
kth Malliavin derivative in a compact form. This is stated as follows: for 
any subset K = {h 1 ,..., h n } of {1, ... , k}, one sets j(K) =j hl ,... ,j h and 
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r(K) = r/ ll , . . . , 77, . Then, one defines 

„(*•'•! a): Di\-±A\{X S ) 



xD^ T 1 \x^---D j{ ^\x k " 

r(h) s r(I„) s 



and 



= Y,d kl ---d kv B\X s ) 



x rPftlx* ■ ■ ■ D j{ l I »\x k s » , 

r(h) s r(I u ) s > 

where in both cases the sum is extended to the set of all partitions of 
{1, . . . , k} = I\ U ■ • ■ U I v . Finally, one sets ot l As) = A l AX s ). Making use of 
this notation, it is shown that the equation satisfied by the fcth derivative 
reads as 

k 

D k:i\ X i = E a J j j ., ,,('••'•!• • • ■ ,r £ -i,r £+1 , ...,r k ) 

e=l 

(A.l) + f (p i Ju ^ Jk {s,r 1 ,...,r k )ds 

J ri V-Vr k 

+ a \ji,...,j k ( s i r ^>---> r k)dW l s ), 

if t > r\ V • ■ ■ V r k , and Dr{l'.'.'.f r k k Xl = otherwise. We prove (2.16) by in- 
duction. The estimate is true for k = 1, with 7i jP = 2Ci iP : this simply fol- 
lows with an application of Burkholder's inequality and Gronwall's lemma 
to (A.l) taken for k = 1. Let us suppose that (2.16) is true up to k — 1. 
As done for k = 1, we apply Burkholder's inequality to (A.l), and, setting 
r = r\ V • ■ • V r k , we get 



n\D^ : -f k xi\p] 

< C k J Y, E Hedl,...,je-l,je + l,-..,3k( r ^ r ^ " " " > ^-1^+1, ■ ■ ■ ,T k )\ P ] 

U=i 

+ (t - rf' 2 ' 1 

rt [ / 

(t-r) 1 / 2 \d kl ---d kv B i (X s )\ 



x £ / E 

hU-Ulu ' 
card(/)<fc-l 



r 



(A.2) 



ds 
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d \P 

+ y £ t \d kl -d hv A\(X a )\) 
1=1 J 

x \D j( / T 1 \x^---D j fr\xH 

I r(h) s r (J„) s I 

+ (t - rf' 2 - 1 J\ [{t - r) l ' 2 \d k B\X s )\ 

d \p 



ds }, 



where, in the last line, we have isolated the term depending on Dl\^'''^ r k k X . 

To estimate the second term in (A.2) we notice that, for any partition 
I\ U • • • U I v of {1, . . . , A;} such that card(I) < k — 1, using (2.16) up to order 
k — 1 we have 



E 



{t-r) 1 ' 2 \d kl ---d kv B\X s )\ 



(A.3) 



+ Y,\dk x ---d K A\{Xs)\ ) \D%^ ■ ■ ■ D j r ^X k ^ 



i=i 



(i) 



<C{|A|t 2 (i 1/2 |5U + I^U) x,t }%(i) Afe - 
where we have defined 



A 



kl := SU P {7card(/i),p H ^7card(/,),p} 

hU-Ul v ={l,...,k} 

card(J)<ifc-l 



and 



Xk = 1 + £(card(J { ) + l) 2 



It is easy to see that 



Xk<{k + lf 



since 



± ca rd (/,) 2 = (* - If ± 2^*2 < <* - 1)> ± ^ = 



fe(fc-l) 5 



- (k-1) 
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so that 

V V 

Xk = l + Y^ card(/z) 2 + 2 ^ card(lj) + v 
1=1 i=i 

<l + (k-l)k + 2k + u 

< 1 + k 2 -k + 3k = {k + lf. 

To estimate the first term in (A. 2), notice that we have as well 

, K . H\^ e j 1 ,...d e - U js + i,.,j k ( r ^ri, ■ ■ .,r e -i,r e+1 , . . . ,r k )\ p ] 

(A.4) 

< C{]A|t_ 1 (ry2] jB | fc _ 1 + \A\ k ^fre p (t) x Z, 

with 

\p := SU P {7card(7i),p H ^7card(^),p}- 

iiU---U/„={l,...,fc-l} 

We remark that a[ and A^ are defined by means of the 7's up to order 
k-1. 

Collecting (A.2), (A.3) and (A.4), we get 

< ^ >P |l^l£i(* 1/2 |5U-i + I^U-i) fc2p e P (t) A ^ 

+ (t- ryl 2 \A\ k k-2{t l/2 \B\k + \A\ k )^ 2 r>e p {t) X ^ 

m r-t \ 

+ (t- r yi^(t i i 2 \B\ l + \A\ l yY J / ®\ D k:zi k k xi\ p ds\ 

k=l Jr J 

< c k ,M\ti(t 1/2 \B\k + I^U) (fc+1)2 %(0 A ^ VAi - 

x (1 + C fcjP t?/ 2 (t 1/2 |£|i + I^Kf e p (tf^) 

< C k>p \A\i P -i(t 1/2 \B\k + |^| fc )( fc+1 ) 2 %(t) A ^ vA ^ +2 ^-, 

where we have applied Gronwall's lemma to get the second inequality. The 
constant C k , P may vary from line to line, but never depends on t nor on 

the bounds on B and A. We recursively define 7^ by setting 7^ := V 

(2) 

X k + 2C\ P , and we finally obtain (2.16). □ 
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Proof of Lemma 2.2. Step 1. We first use the decomposition D s X t = 
Y t Z s A(X s ) (see, e.g., [16]) and write 

a Xt =Y t f Z s A(X s )A(X s )*Z*dsY t * 

(A.5) J ° ( 

= Y t U t Y t *, 

where we have set Ut = J Z S A(X S )A(X S )* Z* ds. Notice that U% is a positive 
operator, and that for any £ G M m we have 

<£,^£> = [ t (A(x 8 yz* a z,A(x 8 yz* a Od S = V f (z s a,(x s ),o 2 . 

Jo j=1 Jo 

From identity (A.5) it follows that det a Xt = (det Y t ) 2 det f/ t = (det Z t )~ 2 det J7 t . 
Hence, applying Holder's inequality, 

E[|detax t n < (E[|detZ t |^]E[(det^)" 2p ]) 1/2 

(A.6) 

<c P u4 P {t) m n{^u t r 2p \) 1/ \ 

where in the last step we have used bound (2.15) on the entries of Zt- 
Step 2. Let X t = inf|£i =1 (£, Ut£) be the smallest eigenvalue of Ut, so that 

E{(detU t )~ 2p } < E[A t " 2mp ]. We evaluate P(A t < e). 

For any £ such that |£| = 1, using the elementary inequality (a + 6) 2 > 
a 2 /2 - b 2 we get 

d 1 d d 

J2(Z s Aj(X s ),0 2 > ^(Aj(x),0 2 - J2(Z s Aj(X s ) - Mx),0 2 
j=l j=l j=l 

1 d 

>-c*-J2\z s M x s)-Mx)\ 2 , 

where in the last step we have used the ellipticity assumption (E). For any 
e > and a > such that as <t, the previous inequality gives 



[ \ac*e- sup ^ ae^2\Z s Aj(X s ) - Aj(x)\ 2 \ < e 
3=1 J / 



s<ae 



and thus, if we take a = 4/c* in order to have ac*/2 = 2 and apply Markov's 
inequality, we obtain 



P(A t < e) < P^sup |^ |Z^(A%) - A,(x)| 2 j > ^ J 

.„ d 
3=1 



(A.7) 
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where the last holds for all q > 1. Now, to estimate the last term we claim 
that, for all j = 1, ... , d, 

(A.8) E\ B uT>\Z a Aj{X a ) - Aj(x)M < Ctf(lV 2 \B\ \A\l + |A| 2 ) 2 «ef ? (i) c , 

l s<t J 

for a constant C depending on q, m, d but not on the bounds on B and A. 
From (A. 7) and this last estimate, it follows that 

P(A* < e) < C w 4(* 1/2 l^lol^li + \A\l) 2 *eZ q (tf^, 

for any e such that 4e/c* < 1 A t. 

Step 3. We finally estimate E[A i T 2mp ]. We write 

oo 

E[A7 2mp ] = n\; 2mp l {Xt>1} ] + ^E[A7 2 ^l {1/(fc+1)<At < 1/fc} ] 

k=l 

oo 

< 1 + J2( k + l) 2mp P(V(fc + 1) < A t < 1/fc), 
fe=l 

and separate the contribution of the sum over k > to obtain 

E[\; 2mp ] < 1 + ^ (Jfe + l) 2mp P(l/(A; + 1) < X t < l/k) 

l<fc<4/(te») 

+ Yl (k + l) 2mp n\t<l/k) 

fc>4/(tc*) 

2mp 



<!+(— + !) P(A t < 1) 



+ ^, m , d ef 9 (t) c — *±-(i}l 2 \B\ Q \A\l + \A\ 2 ) 2 * 



fc>4/(tc») 

We finally take q = 2mp + 2 in order to get convergent series. This last 
estimate, together with (A. 6), gives the desired result. 

PROOF of ( A.8) . We apply Ito's formula to the product Z t Aj (X t ) and get 
d(Z t Aj(X t )) = Zt^idAjB - dBA) + Y^dA l {dA l A J - dAjA^Xjdt 

(A.9) + Zt (±d kl d k2 A j A^A'tA (X t ) dt 

d 

+ Z t J2(8A j A l - dAiAj^X^dW^t). 



i=i 
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Hence, by Burkholder's inequality, 



sup E sup\(Z s Aj(X s ) - Aj(x)y 



i=l,...,m L s<t 



< C^dBlol^li + \B\i\A\o + \A\l\A\ + \A\ 2 \A\ 2 ) 2 i 

+ t i- 1 (\A\ 1 \A\o) 2q } 

x / e[ sup \{Z u ) htk \ 2q \du 



JO 4i,.fe=l,... ) m 



<CV(& 2 \B\ \A\l + \A\l)*4.(t) 



where the constant C depends on q, m and d, but not on t and on the bounds 
on B and A and their derivatives. In the last step, we have once again used 
bound (2.15) on the entries of Z. □ 

A. 2. Proof of Proposition 3.1. We first collect the basic facts we need 
to give the proof of Proposition 3.1. We will start by proving existence and 
uniqueness of strong solutions for the following equation: 



whose coefficients are defined on the whole real line [a,b and 7 are the 
functions appearing in (3.1)]. Once we have established that the unique 
strong solution of (A. 10) is a.s. positive, then (A. 10) will coincide with the 
original equation (3.1). 

The proof of Proposition 3.1 is split in the following two short lemmas. 

Lemma A.l. Assume condition (sO) of Proposition 3.1. Then existence 
and uniqueness of strong solutions hold for (A. 10). Moreover, for any initial 
condition x > the solution is a.s. positive, P(A"j > 0;t > 0) = 1. 

Proof. Existence of nonexplosive weak solutions for (A. 10) follows from 
continuity and sub-linear growth of drift and diffusion coefficients. The ex- 
istence of weak solutions together with pathwise uniqueness imply the exis- 
tence of strong solutions (cf. [12], Proposition 5.3.20 and Corollary 5.3.23). 
Pathwise uniqueness follows in its turn from a well-known theorem of unique- 
ness of Yamada and Watanabe (cf. [12], Proposition 5.2.13). Indeed, as 
a, 6, 7 S C\ the diffusion coefficient of (A. 10) is locally Holder-continuous 
of exponent a > 1/2 and the drift coefficient is locally Lipschitz-continuous. 



(A.10) 
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We apply the standard localization argument for locally Lipschitz coeffi- 
cients and the Yamada-Watanabe theorem to establish that solutions are 
pathwise unique up to their exit time from a compact ball, and hence path- 
wise uniqueness holds for (A. 10). □ 

Lemma A. 2 deals with the second part of Proposition 3.1, that is, the 
behavior at zero. The proof is based on Feller's test for explosions of solutions 
of one-dimensional SDEs (cf. [12], Theorem 5.5.29). Letting r denote the 
exit time from (0, oo), that is, r = inf{t > 0:Xt ^ (0, oo)} with inf = oo, 
we have to verify that 

(A. 11) \\m.p c {x) = — oo 

x— >o 

with p c defined by 

(A.12) Pc (x):=J^ exp (-2 £ a ^ ~ b ^ Z dz) dy, x > 0, 

for a fixed c > 0. Property (A. 11) implies that P(r = oo) = 1, then r = tq 
with To as defined in Proposition 3.1, because the solution of (A. 10) does not 
explode at oo (cf. Lemma A.l). The inner integral in (A.12) is well defined 
and finite for any y > because "f(z) 2 > for any z > and 7 is continuous. 



Remark A.l. The conclusion does not depend on the choice of c € (0, 00) . 

Lemma A. 2. Assume (s0), and let X = (Xt;t > 0) denote the unique 
strong solution of (A. 10) for initial condition x > 0. Then the statements of 
Proposition 3.1 on the stopping time tq hold true. 

Proof. We prove (A. 11), for c = 1. We assume without restriction that 
x < 1 and distinguish the two cases. 

Case a > 1/2. We have a{z) > a(0) - \a\\z, z>0. Then 

a(z) - b(z)z > a(0) - (|a|i + b(z))z 
■-f(z) 2 z 2a ~ -)(z) 2 z 2a 

a(0) |a|i + |&|o 



> 



\j\ 2 z 2a 7 (*) 2 z 2 "- 1 ' 

is integrable at zero by (si)', and then there exists a positive 
constant K such that 

y a(z) - b{z)z 2o(0) f 1 dz „ 

> , ,0 / ~^TZ + & 



1)+K 



7 (Z) 2 Z 2 - - | 7 |§ J y Z 2 « 



Jy 
2a(0) / 1 



(2a-l)| 7 |gVl/ 



,2o;-l 
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hence 

Ji t 2 P \(2a-l)h\ 2 J 



Case a = 1/2. By (s2), 



2a(z) > 1 



7(z) 2 z z' 



for z < x. Hence 



^{z) 2 z z 7(z) 2 ' 
and thus, being integrable at zero, for x <x we have 



Pi{x) < —C J exp^y -dz^jdy 



1 l 

-C / - dy — > — oo. 



□ 
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